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(54) Methods and apparatus for calibrating ct x-ray beam tracking loop 



(57) The present invention is, in one embodiment, a 
method for determining tracking control parameters for 
positioning an x-ray beam (1 6) of a computed tomogra- 
phy imaging system (10) having a movable collimator 
(52) positionable in steps and a detector array (20) 
including a plurality of rows (54, 56, 58, 60) of detector 
elements. The method includes steps of obtaining 
detector samples at a series of collimator step positions 
while determining a position of a focal spot of the xray 
beam (1 6); determining a beam position for each detec- 
tor element at each collimator step utilizing the deter- 
mined focal spot positions, a nominal focal spot length, 
and geometric parameters of the x-ray beam (1 6), colli- 
mator (52) , and detector array; and determining a cali- 
bration parameter utilizing information so obtained. For 
example, in determining a target beam position at which 
to maintain the x-ray beam (1 6), a detector element dif- 
ferential error is determined according to ratios of suc- 
cessive collimator step positions; and a target beam 
position is selected for an isocenter element in accord- 
ance with the determined element differential errors. 
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Description 

[0001] This invention relates generally to computed 
tomography (CT) imaging and, more particularly, to 
methods and apparatus for calibration of z-axis tracking 5 
loops for positioning a CT x-ray beam of a multi-slice CT 
imaging system. 

[0002] In at least one known computed tomography 
(CT) imaging system configuration, an x-ray source 
projects a fan-shaped beam which is collimated to lie 10 
within an X-Y plane of a Cartesian coordinate system 
and generally referred to as the "imaging plane". The x- 
ray beam passes through the object being imaged, such 
as a patient. The beam, after being attenuated by the 
object, impinges upon an array of radiation detectors. 15 
The intensity of the attenuated beam radiation received 
at the detector array is dependent upon the attenuation 
of the x-ray beam by the object. Each detector element 
of the array produces a separate electrical signal that is 
a measurement of the beam attenuation at the detector 20 
location. The attenuation measurements from all the 
detectors are acquired separately to produce a trans- 
mission profile. 

[0003] In known third generation CT systems, the x- 
ray source and the detector array are rotated with a gan- 25 
try within the imaging plane and around the object to be 
imaged so that the angle at which the x-ray beam inter- 
sects the object constantly changes. A group of x-ray 
attenuation measurements, i.e., projection data, from 
the detector array at one gantry angle is referred to as a 30 
"view". A "scan" of the object comprises a set of views 
made at different gantry angles, or view angles, during 
one revolution of the x-ray source and detector. In an 
axial scan, the projection data is processed to construct 
an image that corresponds to a two-dimensional slice 35 
taken through the object. One method for reconstructing 
an image from a set of projection data is referred to in 
the art as the filtered back projection technique. This 
process converts the attenuation measurements from a 
scan into integers called "CT numbers" or "Hounsfield 40 
units", which are used to control the brightness of a cor- 
responding pixel on a cathode ray tube display. 
[0004] In a multi-slice system, movement of an x- 
ray beam penumbra over detector elements having dis- 
similar response functions can cause signal changes 45 
resulting in image artifacts. Opening system collimation 
to keep detector elements in the x-ray beam umbra can 
prevent artifacts but increases patient dosage. Known 
CT systems utilize a closed-loop z-axis tracking system 
to position the x-ray beam relative to a detector array. It so 
would be desirable to provide improved methods and 
apparatus for calibration of such systems. In particular, 
it would be desirable to provide improved methods and 
apparatus for determining calibration parameters such 
as: (1 ) a target beam position at which to maintain the x- 55 
ray beam; (2) a transfer function to convert sensed 
tracking information into a beam position in millimeters; 
and (3) valid limits of the transfer function. 



[0005] There is therefore provided, in one embodi- 
ment, a method for determining tracking control param- 
eters for positioning an x-ray beam of a computed 
tomography imaging system having a movable collima- 
tor positionable in steps and a detector array including a 
plurality of rows of detector elements. The method 
includes steps of obtaining detector samples at a plural- 
ity of collimator step positions while determining a posi- 
tion of a focal spot of the x-ray beam; determining a 
beam position for each detector element at each colli- 
mator step utilizing the determined focal spot positions, 
a nominal focal spot length, and geometric parameters 
of the x-ray beam, collimator, and detector array; and 
determining a calibration parameter utilizing information 
so obtained. For example, in determining a target beam 
position at which to maintain the x-ray beam, the 
method also includes steps of determining an detector 
element differential error according to ratios of succes- 
sive collimator step positions; and selecting a target 
beam position for an isocenter element in accordance 
with the determined element differential errors. 
[0006] The above described system provides 
improved tracking calibration for CT imaging systems 
utilizing z-axis tracking loops for positioning x-ray 
beams. 

[0007] The invention will now be described in 
greater detail by way of example, with reference to the 
drawings in which: 

Figure 1 is a pictorial view of a CT imaging system. 

Figure 2 is a block schematic diagram of the system 
illustrated in Figure 1 . 

Figure 3 is a schematic view of a portion of the CT 
imaging system shown in Figure 1 showing an 
embodiment of a z-axis position system of the 
present invention. 

Figure 4 is a flow diagram an embodiment of a z- 
axis tracking loop of the present invention. 

Figure 5 is a flow diagram of a method for calibrat- 
ing tracking loop parameters. 

[0008] Referring to Figures 1 and 2, a computed 
tomograph (CT) imaging system 10 is shown as includ- 
ing a gantry 1 2 representative of a "third generation" CT 
scanner. Gantry 12 has an x-ray source 14 that projects 
a beam of x-rays 16 toward a detector array 18 on the 
opposite side of gantry 12. Detector array 18 is formed 
by detector elements 20 that together sense the pro- 
jected x-rays that pass through an object 22, for exam- 
ple a medical patient. Each detector element 20 
produces an electrical signal that represents the inten- 
sity of an impinging x-ray beam and hence the attenua- 
tion of the beam as it passes through patient 22. During 
a scan to acquire x-ray projection data, gantry 12 and 
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the components mounted thereon rotate about a center 
of rotation or isocenter 24. 

[0009] Rotation of gantry 12 and the operation of x- 
ray source 14 are governed by a control mechanism 26 
of CT system 10. Control mechanism 26 includes an x- 5 
ray controller 28 that provides power and timing signals 
to x-ray source 14 and a gantry motor controller 30 that 
controls the rotational speed and position of gantry 12. 
A data acquisition system (DAS) 32 in control mecha- 
nism 26 samples analog data from detector elements 10 
20 and converts the data to digital signals for subse- 
quent processing. An image reconstructor 34 receives 
sampled and digitized x-ray data from DAS 32 and per- 
forms high-speed image reconstruction. The recon- 
structed image is applied as an input to a computer 36 15 
that stores the image in a mass storage device 38. 
[0010] Computer 36 also receives commands and 
scanning parameters from an operator via console 40 
that has a keyboard. An associated cathode ray tube 
display 42 allows the operator to observe the recon- 20 
structed image and other data from computer 36. The 
operator supplied commands and parameters are used 
by computer 36 to provide control signals and informa- 
tion to DAS 32, x-ray controller 28 and gantry motor 
controller 30. In addition, computer 36 operates a table 25 
motor controller 44 that controls a motorized table 46 to 
position patient 22 in gantry 12. Particularly, table 46 
moves portions of patient 22 through gantry opening 48. 
[0011] In one embodiment, and as shown in Figure 
3, x-ray beam 1 6 emanates from a focal spot 50 of x-ray 30 
source 14 (Figure 2). X-ray beam 16 is collimated by 
collimator 52, and collimated beam 16 is projected 
toward detector array 1 8. Detector array 1 8 is fabricated 
in a multi-slice configuration and includes detector ele- 
ment rows 54, 56, 58 and 60 for projection data collec- 35 
tion. A plane 86, generally referred to as the "fan beam 
plane", contains the centerline of focal spot 50 and the 
centerline of beam 16. Fan beam plane 86 is illustrated 
in Figure 3 as being aligned with a centerline D 0 of 
detector array 18, although fan beam plane 86 will not 40 
always be so aligned. Detector element rows 62, 64, 66 
and 68 serve as z-position detectors for determining a 
z-axis position of x-ray beam 16. In one embodiment, 
detector rows 62, 64, 66, and 68 are rows of detector 
array 18. Outer rows 62 and 68 are selected to be at 45 
least substantially within penumbra 70 of beam 16. 
Inner rows 64 and 66 are selected to be at least sub- 
stantially within umbra 72 of beam 1 6. "At least substan- 
tially within" means either entirely within or at least 
sufficiently within so that outer row 62 and 68 signal 50 
intensities depend on an x-ray beam position and inner 
row 64 and 66 signal intensities provide references 
against which outer row signals are compared. In one 
embodiment, collimator 52 includes tapered cams 74 
and 76. (Where it is stated herein that a cam "has a 55 
taper," it is not intended to exclude cams having a taper 
of zero unless otherwise stated.) X-ray controller 28 
controls positioning of cams 74 and 76. Each cam can 



be independently positioned to alter position and width 
of x-ray umbra 72 relative to an edge (not shown) of 
detector array 1 8. 

[0012] As shown in Figure 4, one embodiment of a 
closed-loop method for positioning beam 16 comprises 
comparing signals representative of x-ray intensity 
received from different rows of detector elements and 
positioning an x-ray beam in accordance with results of 
the comparison. In one embodiment, signals represent- 
ative of x-ray intensity from detector rows 62, 64, 66 and 
68 are summed 78 to obtain row sums. The summation 
is over views taken in a 20-millisecond interval. For 
example, after the analog signals are converted to dig- 
ital format, hardware circuitry (not shown) in DAS 32 
performs offset correction and determines row sums 
from signals received from outer row 62 and from inner 
row 64. A corrected ratio R is determined 80 by deter- 
mining a ratio of a sum of signals received from outer 
row 62 to a sum of signals received from inner row 64 
and multiplying the ratio by a ratio correction factor. The 
ratio correction factor, determined from imaging system 
1 0 calibration, accounts for different relative DAS gains 
between outer row 62 and inner row 64. 
[0013] Beam position Z(R) then is determined 82, 
in millimeters relative to a centerline. Beam position Zis 
obtained by applying a predetermined beam position 
transfer function to the corrected ratio to calculate the x- 
ray beam position. The beam position transfer function 
Z(R) is represented, for example, by a fourth-degree 
polynomial having predetermined coefficients: 

Z(R) = a + bR + cR 2 + dR 3 + eR 4 

[0014] Beam position transfer function Z(R) and its 
limits are specified at imaging system 10 calibration. 
[0015] A new collimator position is then determined 
84. A focal spot position f is determined 84 from beam 
position Z, current collimator position C and other sys- 
tem 10 geometric parameters in accordance with: 

(Z-C-7 Z ) 

f = -1 71 - T + C + T 7 

fm zz (f fs ) 

where T z represents a current taper of cam 74, fm zz 
represents a focal spot magnification factor at rows 62 
and 64 and is a function of focal spot size, and l fs repre- 
sents focal spot 50 length. A new position for collimator 
52 then is determined 84 for a detector element 20 posi- 
tioned toward isocenter 24. Collimator 52 is reposi- 
tioned where an edge (not shown) of collimator 52 
would meet a line between focal spot position f and a 
target beam position Z f which has been specified at 
imaging system 10 calibration. New collimator position 
C n thus is determined in accordance with: 
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(Z t -f) tem 10 to more quickly detect an end of the blockage 

cm y — j + f and to resume closed-loop positioning. 

[0019] Expected signal Sx is written as: 



where cm } represents a current collimator magnification 5 
factor at detector element 20 positioned toward iso- 
center 24 and is a function of focal spot size, and l fs rep- 
resents focal spot 50 length. 

[0016] In one embodiment, steps 78, 80, 82, and 84 
are performed independently for each side of collimator 10 
52 at intervals to continuously obtain new positions for 
each side of collimator 52. These intervals are, in one 
embodiment, 20 milliseconds, to sample the x-ray beam 
16 position 25 times during a 0.5 second scan to mini- 
mize control loop lag error. However, in other embodi- 15 
ments, the interval is between 5 milliseconds and 50 
milliseconds. In still other embodiments, the interval is 
between a minimum value sufficient to avoid effects of 
quantum noise and high frequency variation (such as 
due to x-ray tube anode movement at a run frequency 20 
between 50 Hz and 160 Hz) and a maximum contrained 
by a slew rate of the sag curve. Sampling the changing 
sag curve frequently avoids excessive positioning error. 
(Sag is a periodic movement of x-ray beam 1 6 resulting 
from gravity and from centrifugal forces acting on 25 
mechanical structure during a rotation of gantry 12.) 
[0017] During patient scanning, z-position detec- 
tors 62, 64, 66 or 68 may become blocked by patient 
clothing, blankets, or other object. After blockage of a z- 
position detector 62, 64, 66, or 68 has been detected, or 30 
when x-ray source 14 first turns on, the loop sample 
interval is adjusted downward. In one embodiment, the 
loop sample interval is adjusted downward to 5 millisec- 
onds. After 4 milliseconds of stabilization, the position of 
the beam is measured and collimator positioning is 35 
started to further minimize initial position errors. 
[0018] During a blockage, loop operation is sus- 
pended. To determine if any z-position detectors are 
blocked, a signal from a last data detector element 90 
adjacent a z-position detector 62, 64, 66 or 68 is com- 40 
pared to an expected signal Sx. Z-position detector 
blockage is assumed, in one embodiment, if a last data 
detector element 20 signal is less than 0.9 times 
expected signal Sx. In other embodiments, detector 
blockage is assumed when a last data detector element 45 
20 signal is less than a value between 0.95 and 0.5 
times expected signal Sx. (It is desirable to make this 
value as large in magnitude as possible to identify 
patient blockage as quickly as possible, thereby avoid- 
ing mispositioning of x-ray beam 1 6 due to corrupted Z- so 
measurement data. A maximum of 0.95 is used in one 
embodiment because it is known that x-ray scatter 
blockage from large patients 22, for example, can 
reduce a signal to 0.95 times the expected value.) Dur- 
ing a blockage, collimator positioning is suspended. 55 
However, position measurement continues at an inter- 
val that is decreased from 20 to 5 milliseconds. The 
decreased measurement interval allows imaging sys- 



Sx = gmA * csf * t* g, 

where gmA is a generator current mA signal propor- 
tional to an x-ray source 14 energizing current, csf is a 
scale factor determined at system 10 calibration, t is a 
DAS sample time period, and g is a gain factor. Gain 
factor g allows expected signal Sx to be adjusted 
according to a gain value used for scanning. In one 
embodiment, this gain value is selectable from a plural- 
ity of gain values available in system 10. 
[0020] In one embodiment, closed loop tracking is 
suspended when signal corruption is detected. Signal 
corruption is detected, for example, by determining an 
actual focal spot length from a beam position and a col- 
limator position, and comparing the actual focal spot 
length to a nominal focal spot length. When a difference 
of, for example, more than 0.1 millimeter is detected 
between the actual focal spot length and the nominal 
spot length, corruption is assumed to exist and collima- 
tor positioning is suspended. (In other embodiments, a 
difference threshold for assuming corruption is as small 
as 0.05 millimeter or as large as about 0.6 millimeter. In 
still other embodiments, a value is selected between a 
lower limit set by higher probabilities of false activation 
due to noise, x-ray scatter and/or momentary beam 
position disturbances and an upper limit that still pro- 
vides some of the advantages of tracking.) However, 
beam position measurement continues at a decreased 
interval, as when a blockage is detected. Such corrup- 
tion may occur, for example, for a short time just prior to 
or just following detection of a patient blockage. If the 
corruption persists, for example, over 90° of rotation of 
gantry 12 without detecting a patient blockage, a mal- 
function of the tracking system requiring servicing has 
likely occurred. In such an event, a scan is immediately 
aborted to avoid patient dose and collection of non- 
diagnostic quality images. In other embodiments, a limit 
is set from as little as 45° to as much as 360° of a rota- 
tion of gantry 12. In other embodiments, a limit is set 
between a value at which a false alarm rate due to scat- 
ter and/or an occasional exceptionally long partial 
patient 22 blockage is acceptable and an upper limit 
representing a design choice as to how long compro- 
mised operation (high dose and/or non-diagnostic qual- 
ity images) can be tolerated before terminating a scan. 
[0021] After system 10 has been switched off, posi- 
tion of focal spot 50 changes as source 14 cools over 
time. In one embodiment, before system 10 is switched 
on again, an initial focal spot position is approximated 
from information obtained when a focal spot position 
was last measured. An approximation of a linear func- 
tion is used to model focal spot position change during 
cooling in one embodiment, and in another embodi- 
ment, the linear function is a 97 nanometer per second 
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linear function. Because position change with cooling is 
an exponential function, the linear approximation is 
clamped at 0.1 5 millimeters. This clamping corresponds 
to approximately 20% of a cooling change in system 10 
when fully cold, where a linear approximation to the 5 
exponential function suffices. A fully cold position 
requires 8 to 12 hours without patient scanning, and a 
tube warm up prior to patient scanning is normally 
requested if the tube has been off more than 1 hour. 
Therefore, a fully cold position, although possible, is not w 
likely during normal patient scanning. During tube warm 
up a current measured position of the focal spot is 
established again for initial positioning of the collimator. 
[0022] Several tracking loop parameters described 
herein, specifically, beam position transfer function Z(R) 15 
and its limits and target beam position Z t , are deter- 
mined at system 10 calibration. Figure 5 illustrates one 
embodiment of a method for calibrating tracking loop 
parameters. In this embodiment, data from a stationary 
sweep scan is collected 100 while collimator 52 is 20 
stepped through a sequence of z-axis positions. Beam 
16 is incremented 0.3 millimeters on detector array 18 
exposure surface for each collimator 52 step position. 
The sweep scan data is offset-corrected and view aver- 
aged 1 02 to obtain a set of detector samples for each 25 
collimator 52 step position. A position of the focal spot is 
then determined 1 04. A collimator 52 z-axis position off- 
set from detector array centerline D 0 is determined 1 04, 
as the point where outer rows 62 and 68 receive signals 
of half-maximum intensity at full detector element 20 30 
width. Position of focal spot 50 during sweep scan then 
is determined 104 from collimator 52 z-axis offset and 
nominal system 10 geometric parameters. 
[0023] A beam 16 position is determined 106 for 
each detector element 20 at each collimator 52 step 35 
position. Beam 1 6 positions are determined from sweep 
scan focal spot 50 position, nominal length of focal spot 
50, and nominal system 10 geometry. 
[0024] Target beam position Z t then is determined 
1 08 for detector element 20 positioned toward isocenter 40 
24. When beam 16 is directed at target beam position 
Z t , beam 16 is sufficiently close to detector array 18 
edge 92 to prevent imaging artifacts but is far enough 
away to minimize patient dosage. To determine target 
beam position Z t , ratios of detector samples for succes- 45 
sive collimator 52 step positions are utilized to deter- 
mine a detector differential error. A reconstruction error 
sensitivity function w(i) then is applied to weight the 
detector differential error. Reconstruction error sensitiv- 
ity function w(i) is related to the percent positive contri- so 
bution of a detector element 20 as a function of its radial 
distance from isocenter 24. Function w(i), in one 
embodiment, is computed from nominal system geome- 
try. In another embodiment, w{i) is empirically deter- 
mined. For example, the following equations describe 55 
an empirical determination of w(i): 

b(i) = .018, 0 < / < 5 



b(i) = .035 + .00075x(i — 5), 5 < / < 213 
b(i) = .414 + .00365x(i — 21 3), 21 4 < i < n 

where / represents detector element position from iso- 
center 24 and b(i) represents an artifact threshold, i.e. a 
percent differential error, for a double detector element 
20 error. Reconstruction error sensitivity function w(i) 
then is determined in accordance with: 

w(i) = 0.18/ b(i). 

[0025] A collimator 52 step position SP is deter- 
mined for which the weighted detector differential error 
exceeds a limit L empirically known to produce image 
artifacts, for example, 0.04 percent. Target beam posi- 
tion Z t then is set for the isocenter detector element at a 
distance just preceding SP by an amount exceeding 
applicable tracking loop positioning error. 
[0026] Beam position transfer function Z(R) then is 
determined 1 1 0 for a ratio R of an average of outer row 
62 to inner row 64 signals for a set of detector elements 
at an extreme end of x-ray fan beam 1 6. Beam 1 6 posi- 
tions, determined 106 for each collimator 52 step posi- 
tion, are fitted to the ratio for each collimator 52 step 
position with a fourth-degree polynomial, for example, in 
accordance with: 

Z(R) = a + bR + cR 2 + dR 3 + eR 4 

over a suitable ratio range between a maximum and 
minimum for the sequence of steps. 
[0027] A valid position measurement range for Z(R) 
is determined 112 as between end limits of the set of 
collimator 52 step positions for which an error between 
a beam 16 position determined by Z and an actual 
beam 16 position is less than a predetermined limit, for 
example, 0.2 millimeters. In other embodiments, the 
predetermined limit is between 0.1 millimeters to 0.6 
millimeters. In still other embodiments, the predeter- 
mined limit is set at a value between a lower limit just 
above a value at which a range of beam 1 6 position that 
can be precisely measured is too limited, and just below 
a lower limit that is deemed to create tracking errors so 
large as to unacceptably compromise the benefits of 
tracking. 

[0028] The above-described tracking loop senses 
the signal ratio between detector rows and moves sys- 
tem collimation to maintain the x-ray beam very close to 
the imaging system detector array edge during patient 
scanning. As a result, patient x-ray dosage is reduced 
20 to 40 percent without sacrificing image quality. 
[0029] Other functions can be utilized in place of 
beam position transfer function Z(R) and also in place 
of reconstruction error sensitivity function w(i). 
[0030] In some embodiments, the methods 
described herein are implemented by software, 
firmware, or by a combination thereof controlling either 
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computer 36, image reconstructor 34, or both. Also, 
additional z-detector rows can be provided. In such an 
embodiment, various combinations of z-detector row 
signals can be used as the inner and outer row signals, 
thereby becoming identified as such, or a different 5 
and/or more elaborate transfer function can be used to 
determine a beam position. 

[0031] The above described calibration methods 
and apparatus provides improved calibration for z-axis 
tracking loops for positioning x-ray beams on multi-slice 10 
detectors of CT imaging systems. The methods and 
apparatus provide a target beam position at which to 
maintain the x-ray beam, a transfer function to convert 
detector ratio information into a beam position in millim- 
eters (or other suitable units via conversion factors), and 15 
valid limits of the ratio to beam position transfer function. 
[0032] It should be understood that system 10 is 
described herein by way of example only, and the inven- 
tion can be practiced in connection with other types of 
imaging systems. Furthermore, it will be recognized by 20 
those skilled in the art that the calibration system 
described herein is also useful for other applications 
which require x-ray beam tracking calibration, such as 
for object location or sensing of movement. 
[0033] For the sake of good order, various features 25 
of the invention are set out in the following clauses:- 

1 . A method for determining tracking control param- 
eters for positioning an x-ray beam of a computed 
tomography imaging system, the imaging system 30 
including a movable collimator positionable in steps 
and a detector array including a plurality of rows of 
detector elements, said method comprising the 
steps of: 

35 

obtaining detector samples at a plurality of col- 
limator step positions while determining a posi- 
tion of a focal spot of the x-ray beam; 

determining a beam position for each detector 40 
element at each collimator step utilizing the 
determined focal spot positions, a nominal 
focal spot length, and geometric parameters of 
the x-ray beam, collimator, and detector array; 

45 

determining an detector element differential 
error according to ratios of successive collima- 
tor step positions; and 

selecting a target beam position for an iso- 50 
center element in accordance with the deter- 
mined element differential errors. 

2. A method in accordance with Clause 1 wherein 

the plurality of detector rows are z-axis detector 55 
rows, and the detector array has a centerline per- 
pendicular to the z-axis, an outer detector row, and 
an inner detector row; said determining a position 



of a focal spot of the x-ray beam comprises the 
steps of: 

determining a collimator z-axis position offset 
from the detector array centerline at a point at 
which outer detector row signals are reduced to 
a full width at a half maximum; and 

determining a focal spot position as a function 
of the determined collimator z-axis position and 
the geometric parameters of the x-ray beam, 
collimator, and detector array. 

3. A method in accordance with Clause 1 further 
comprising the step of offset-correcting and view- 
averaging the obtained detector samples at a plu- 
rality of collimator step positions to obtain a set of 
detector samples for each collimator step position 
used in said steps of determining a beam position 
transfer function and determining a differential error 
for selection of the target beam position. 

4. A method in accordance with Clause 1 wherein 
selecting a target beam position for an isocenter 
detector element in accordance with the deter- 
mined element differential errors comprises the 
steps of: 

weighting the detector element differential error 
by a reconstruction error sensitivity function; 
determining a step position at which the 
weighted detector element differential error 
exceeds a predetermined limit; and 
setting a tracking beam position for the iso- 
center detector element at a distance from the 
determined step position preceding a step that 
exceeds a predetermined artifact limit by an 
amount that exceeds a tracking loop position- 
ing error. 

5. A method in accordance with Clause 4 wherein 
the reconstruction error sensitivity function is detec- 
tor element dependent. 

6. A method in accordance with Clause 5 wherein 
the reconstruction error sensitivity varies according 
to a distance of the detector element from an iso- 
center element. 

7. A method in accordance with Clause 4 wherein 
the detector rows have at least 214 elements on 
each side of an isochannel element, and the recon- 
struction error sensitivity function is: 

w(i) = 0.18/ b(i)\ 

where: 

/ = a detector element position from an iso- 
center detector element; 
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b(i) = an artifact threshold (% differential error) 
for a-double detector element error; and 

b(i) =.018, 0<-/<5 

5 

b(i) = .035 +.00075x(/ - 5) , 5 <-/ < 213 
b(i) = .41 4 + .00365x(/ - 21 3) , 21 4 <-/ < n. 

8. A method for determining tracking control param- 10 
eters for positioning an x-ray fan beam of a com- 
puted tomography imaging system, the imaging 
system including a movable collimator positionable 

in steps and a detector array including a plurality of 
rows of detector elements including at least an 15 
inner row and an outer row, said method comprising 
the steps of: 

obtaining detector samples at a plurality of col- 
limator step positions while determining a posi- 20 
tion of a focal spot of the x-ray fan beam; 

determining a beam position for each detector 
element at each collimator step utilizing the 
determined focal spot positions, a nominal 25 
focal spot length, and geometric parameters of 
the x-ray fan beam, collimator, and detector 
array; and 

determining a beam position transfer function 30 
for a ratio of an average of detector outer row 
signals to detector inner row signals for a set of 
detector elements at an extreme end of the x- 
ray fan beam in accordance with a selected 
approximation over a selected ratio range 35 
between a minimum and a maximum ratio for 
the plurality of collimator step positions. 

9. A method in accordance with Clause 8 wherein 
determining a beam position transfer function com- 40 
prises the steps of fitting, to a polynomial function, 

the determined beam positions at each step as a 
function of the ratio of an average of detector outer 
row signals to detector inner row signals. 

10. A method in accordance with Clause 9 wherein 45 
the polynomial function is a fourth degree polyno- 
mial. 

1 1. A method in accordance with Clause 8 further 
comprising the step of determining a valid meas- 
urement range of the transfer function as end limits 50 
of the plurality of collimator step positions for which 

an error between beam positions computed using 
the transfer function and an actual beam position is 
less than a predetermined limit. 

12. A method in accordance with Clause 11 55 
wherein the predetermined limit is between 0.1 mil- 
limeters and 0.6 mm. 

13. A method in accordance with Clause 11 



wherein the predetermined limit is 0.2 millimeters. 

14. A computed tomography imaging system com- 
prising an x-ray source, a detector array including a 
plurality of rows of detector elements, and a mova- 
ble collimator positionable in steps and configured 
to collimate and position an x-ray beam produced 
by said x-ray source on said detector array, said 
system configured to: 

obtain detector samples at a plurality of colli- 
mator step positions while determining a posi- 
tion of a focal spot of the x-ray beam; 

determine a beam position for each detector 
element at each collimator step utilizing the 
determined focal spot positions, a nominal 
focal spot length, and geometric parameters of 
the x-ray beam, collimator, and detector array; 

determine an detector element differential error 
according to ratios of successive collimator 
step positions; and 

select a target beam position for an isocenter 
element in accordance with the determined 
element differential errors. 

1 5. A system in accordance with Clause 1 4 wherein 
the plurality of detector rows are z-axis detector 
rows, and the detector array has a centerline per- 
pendicular to the z-axis, an outer detector row, and 
an inner detector row; and said system being con- 
figured to determine a position of a focal spot of the 
x-ray beam comprises said system being config- 
ured to: 

determine a collimator z-axis position offset 
from the detector array centerline at a point at 
which outer detector row signals are reduced to 
a full width at a half maximum; and 

determine a focal spot position as a function of 
the determined collimator z-axis position and 
the geometric parameters of the x-ray beam, 
collimator, and detector array. 

16. A system in accordance with Clause 14 further 
configured to offset-correct and view-average the 
obtained detector samples at a plurality of collima- 
tor step positions to obtain a set of detector sam- 
ples for each collimator step position used in 
determining said target beam position transfer func- 
tion and in determining said differential error for 
selection of said target beam position. 

1 7. A system in accordance with Clause 1 4 wherein 
said system being configured to select a target 
beam position for an isocenter detector element in 
accordance with the determined element differen- 
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tial errors comprises said system being configured 
to: 

weight the detector element differential error by 
a reconstruction error sensitivity function; 5 

determine a step position at which the 
weighted detector element differential error 
exceeds a predetermined limit; and 

10 

set a tracking beam position for the isocenter 
detector element at a distance from the deter- 
mined step position preceding a step that 
exceeds a predetermined artifact limit by an 
amount that exceeds a tracking loop position- 15 
ing error. 

1 8. A system in. accordance with Clause 1 7 wherein 
the reconstruction error sensitivity function is detec- 
tor element dependent. 20 

1 9. A system in accordance with Clause 1 8 wherein 
the reconstruction error sensitivity varies according 
to a distance of the detector element from an iso- 
center element. 

20. A system in accordance with Clause 1 7 wherein 25 
the detector rows have at least 214 elements on 
each side of an isochannel element, and the recon- 
struction error sensitivity function is: 

w(i) =QABI b(i)\ 30 

where: 

/ = a detector element position from an iso- 
center detector element; 35 

b(i) = an artifact threshold (% differential error) 
for a-double detector element error; and 

b(i) =.018, 0<-/<5 40 

b(i) = .035 +.00075x(/ - 5) , 5 <-/ < 213 

b(i) = .414 + .00365x(/ - 213) , 214 <-/ < n. 

45 

21. A computed tomography imaging system com- 
prising an x-ray source, a detector array including a 
plurality of rows of detector elements, and a mova- 
ble collimator positionable in steps and configured 

to collimate and position an x-ray beam produced so 
by said x-ray source on said detector array, said 
system configured to: 

obtain detector samples at a plurality of colli- 
mator step positions while determining a posi- 55 
tion of a focal spot of the x-ray fan beam; 

determine a beam position for each detector 



element at each collimator step utilizing the 
determined focal spot positions, a nominal 
focal spot length, and geometric parameters of 
the x-ray fan beam, collimator, and detector 
array; and 

determine a beam position transfer function for 
a ratio of an average of detector outer row sig- 
nals to detector inner row signals for a set of 
detector elements at an extreme end of the x- 
ray fan beam in accordance with a selected 
approximation over a selected ratio range 
between a minimum and a maximum ratio for 
the plurality of collimator step positions. 

22. A system in accordance with Clause 21 wherein 
said system being configured to determine a beam 
position transfer function comprises the steps of fit- 
ting, to a polynomial function, the determined beam 
positions at each step as a function of the ratio of an 
average of detector outer row signals to detector 
inner row signals 

23. A system in accordance with Clause 22 wherein 
the polynomial function is a fourth degree polyno- 
mial. 

24. A system in accordance with Clause 21 further 
configured to determine a valid measurement 
range of the transfer function as end limits of the 
plurality of collimator step positions for which an 
error between bean positions computed using the 
transfer function and an actual beam position is 
less than a predetermined limit. 

25. A system in accordance with Clause 24 in which 
the predetermined limit is between 0.1 and 0.6 mil- 
limeters. 

26. A system in accordance with Clause 24 in which 
the predetermined limit is 0.2 millimeters. 

Claims 

1. A method for determining tracking control parame- 
ters for positioning an x-ray beam of a computed 
tomography imaging system, the imaging system 
including a movable collimator positionable in steps 
and a detector array including a plurality of rows of 
detector elements, said method comprising the 
steps of: 

obtaining detector samples at a plurality of col- 
limator step positions while determining a posi- 
tion of a focal spot of the x-ray beam; 

determining a beam position for each detector 
element at each collimator step utilizing the 
determined focal spot positions, a nominal 
focal spot length, and geometric parameters of 
the x-ray beam, collimator, and detector array; 
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determining an detector element differential 
error according to ratios of successive collima- 
tor step positions; and 

selecting a target beam position for an iso- 5 
center element in accordance with the deter- 
mined element differential errors. 

2. A method in accordance with Claim 1 wherein the 
plurality of detector rows are z-axis detector rows, 10 
and the detector array has a centerline perpendicu- 
lar to the z-axis, an outer detector row, and an inner 
detector row; said determining a position of a focal 
spot of the x-ray beam comprises the steps of: 

15 

determining a collimator z-axis position offset 
from the detector array centerline at a point at 
which outer detector row signals are reduced to 
a full width at a half maximum; and 
determining a focal spot position as a function 20 
of the determined collimator z-axis position and 
the geometric parameters of the x-ray beam, 
collimator, and detector array. 

3. A method in accordance with Claim 1 or 2 further 25 
comprising the step of offset-correcting and view- 
averaging the obtained detector samples at a plu- 
rality of collimator step positions to obtain a set of 
detector samples for each collimator step position 
used in said steps of determining a beam position 30 
transfer function and determining a differential error 

for selection of the target beam position. 

4. A method for determining tracking control parame- 
ters for positioning an x-ray fan beam of a com- 35 
puted tomography imaging system, the imaging 
system including a movable collimator positionable 

in steps and a detector array including a plurality of 
rows of detector elements including at least an 
inner row and an outer row, said method comprising 40 
the steps of: 

obtaining detector samples at a plurality of col- 
limator step positions while determining a posi- 
tion of a focal spot of the x-ray fan beam; 45 

determining a beam position for each detector 
element at each collimator step utilizing the 
determined focal spot positions, a nominal 
focal spot length, and geometric parameters of 50 
the x-ray fan beam, collimator, and detector 
array; and 

determining a beam position transfer function 
for a ratio of an average of detector outer row 55 
signals to detector inner row signals for a set of 
detector elements at an extreme end of the x- 
ray fan beam in accordance with a selected 



approximation over a selected ratio range 
between a minimum and a maximum ratio for 
the plurality of collimator step positions. 

5. A method in accordance with Claim 4 wherein 
determining a beam position transfer function com- 
prises the steps of fitting, to a polynomial function, 
the determined beam positions at each step as a 
function of the ratio of an average of detector outer 
row signals to detector inner row signals. 

6. A method in accordance with Claim 5 wherein the 
polynomial function is a fourth degree polynomial. 

7. A computed tomography imaging system compris- 
ing an x-ray source, a detector array including a plu- 
rality of rows of detector elements, and a movable 
collimator positionable in steps and configured to 
collimate and position an x-ray beam produced by 
said x-ray source on said detector array, said sys- 
tem configured to: 

obtain detector samples at a plurality of colli- 
mator step positions while determining a posi- 
tion of a focal spot of the x-ray beam; 

determine a beam position for each detector 
element at each collimator step utilizing the 
determined focal spot positions, a nominal 
focal spot length, and geometric parameters of 
the x-ray beam, collimator, and detector array; 

determine an detector element differential error 
according to ratios of successive collimator 
step positions; and 

select a target beam position for an isocenter 
element in accordance with the determined 
element differential errors. 

8. A system in accordance with Claim 7 wherein the 
plurality of detector rows are z-axis detector rows, 
and the detector array has a centerline perpendicu- 
lar to the z-axis, an outer detector row, and an inner 
detector row; and said system being configured to 
determine a position of a focal spot of the x-ray 
beam comprises said system being configured to: 

determine a collimator z-axis position offset 
from the detector array centerline at a point at 
which outer detector row signals are reduced to 
a full width at a half maximum; and 

determine a focal spot position as a function of 
the determined collimator z-axis position and 
the geometric parameters of the x-ray beam, 
collimator, and detector array. 



9 



17 EP 1 094 315 A1 18 



9. A computed tomography imaging system compris- 
ing an x-ray source, a detector array including a plu- 
rality of rows of detector elements, and a movable 
collimator positionable in steps and configured to 
collimate and position an x-ray beam produced by 5 
said x-ray source on said detector array, said sys- 
tem configured to: 

obtain detector samples at a plurality of colli- 
mator step positions while determining a posi- 10 
tion of a focal spot of the x-ray fan beam; 

determine a beam position for each detector 
element at each collimator step utilizing the 
determined focal spot positions, a nominal 15 
focal spot length, and geometric parameters of 
the x-ray fan beam, collimator, and detector 
array; and 

determine a beam position transfer function for 20 
a ratio of an average of detector outer row sig- 
nals to detector inner row signals for a set of 
detector elements at an extreme end of the x- 
ray fan beam in accordance with a selected 
approximation over a selected ratio range 25 
between a minimum and a maximum ratio for 
the plurality of collimator step positions. 

10. A system in accordance with Claim 9 wherein said 
system being configured to determine a beam posi- 30 
tion transfer function comprises the steps of fitting, 

to a polynomial function, the determined beam 
positions at each step as a function of the ratio of an 
average of detector outer row signals to detector 
inner row signals 35 
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FIG. 3 
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FIG. 4 



1 . For each row sum signals in all Z channel rows C 
and sum signals over a 20 ms interval of views J 



78 



2. Compute the outer to inner row ratio and Multiply the ratio 
by a DAS gain dependent ratio correction factor 
(from calibration) 



R= 



Outer _row 
Inner _rovo 



ratio _correct ion 



3. Compute beam position by expanding the ratio 
with a 4th degree polynomial (from calibration) 

Z=a+bR-hcR 2 -hdR 3 +eR 4 



J 



80 



82 



84 

4. Using system geometric parameters compute the new collimator position , 

(ZT-f) V 
X C ™<1(I».) J 
* (Z- C- T 7 ) 

Zt = target position 
Cm = collimator mag factor at isocenter, function of fs size 
/ = focal spot position 
C/ = current cam position (at isocenter) 
T z = cam taper (between isocenter and Z channel) 
f 171 xx = focal spot mag factor at Z channel.function of fs size 
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FIG. 5 j r ioo 

Collect a stationary x-ray scan while stepping a movable collimator through a 
sequence of Z-axis locations where the beam is incremented about 0.3 mm on 
the detector for each step 



Offset correct and view average data from the stationary x-ray scan to obtain 
detector samples for each collimator step position 



^-104 



Determine a position of a focal spot during the stationary x-ray scan by 
computing a collimator Z-axis position offset from a detector centerline at a 
point at which outer detector row signals are reduced to full width at half 
maximum. Determine focal spot position from collimator offset thus 
determined using nominal system geometric parameters 

| , -106" 



Determine a beam position for each detector element at each collimator step 
location from the measured focal spot position utilizing a nominal focal spot 
length and nominal system geometry 

\ ^108 
Determine a target beam position Z by: 

• determining a differential error produced by a ratio of signals from successive step 
positions 

• weighting the differential error by a reconstruction error sensitivity function 

• determining a step position at which the weighted differential error goes above limits 
known to produced image artifacts (e.g., 0.04%) 

• setting a tracking target beam position for an isocenter detector element at a 
distance from the step position just preceding the step that exceeds the artifact limit 
by an amount that exceeds any tracking loop positioning error 

1 ~Em 

Determine a beam position transfer function Z (R) for a ratio of an average of outer row 
signals to an average of inner row signals for a set of detector elements at extreme 
ends of the x-ray fan beam by fitting the beam position determined for each step to a 
ratio at each step (e.g., with a 4th degree polynomial) over a suitable ratio range 
between a maximum and a minimum for a sequence of steps 



Determine a valid position measurement range of the transfer function as end limits of 
the set of step positions, where an error between a beam position produced by the 
transfer function and an actual beam position is less than a predetermined limit (e.g., 
0.2 mm) 



14 



EP 1 094 315 A1 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 00307005.9 



Category 



A 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Citation of document with indication, where appropriate, 
of relevant passages 



US 5579359 A 
(TOTH) 

26 November 1996, 

figs. 3-8, column 2, 
line 60 - column 3, line 16, 
column 5, line 7 - column 6, 
line 46, claims 1-7 



DE 19650528 Al 

(GENERAL ELECTRIC CO) 

26 June 1997, 

figs. 1-5, column 2, 
line 28 - column 3, 
column 3, line 26 - 
line 17. 



line 4, 
column 8 , 



WO 00/30538 Al 
(GENERAL ELECTRIC COMPANY) 
02 June 2000, 

fig. 6, page 8, paragraph 2 

page 12, line 2, 

claims 1-4,9,10. 



The present search report has been drawn up for all claims 



Relevant 
to claim 



1,7 



4,5,9, 

10 

4,5,9, 

10 



1,2,7, 

8 

1,2,4, 
7,9,10 



CLASSIFICATION OF THE 
APPLICATION (Int. CI. 7) 



A61B6/03 
G01N23/04 



TECHNICAL FIELDS 
SEARCHED (Int. CI. 7) 



A61B 
G01N 
G21K 
H05G 



Place of search 

VIENNA 



Date of completion of the search 

21-12-2000 



Examiner 

WENNINGER 



CATEGORY OF CITED DOCUMENTS 
X : particularly relevant if taken alone 
Y : particularly relevant if combined with another 

document of the same category 
A : technological background 
O : non-written disclosure 
P : intermediate document 
T : theory or principle underlying the invention 



E : earlier patent document, but published on, or 

after the filing date 
D : document cited in the application 
L : document cited for other reasons 



member of the same patent family, corresponding 
document 



15 



EP 1 094 315 A1 



ANNEX TO THE EUROPEAN SEARCH REPORT 
ON EUROPEAN PATENT APPLICATION NO. EP 00307005.9 



This annex lists the patent family members relating to the patent documents cited in the above-mentioned search report. 
The members are as contained in the EPIDOS INPADOC file on 04.01.2001. The European Patent Office is in no way 
liable for these particulars which are merely given for the purpose of information. 



Patent document cited 


Publication 




Patent family 


Publication 


in search report 


date 




member(s) 


date 


US A 5579359 


26-11-1996 


DE 


Al 


19650527 


26-06-1997 






IL 


AO 


119775 


18-03-1997 






IL 


Al 


119775 


22-09-1999 






JP 


A2 


9285462 


04-11-1997 


DE Al 19650528 


26-06-1997 


IL 


AO 


119772 


18-03-1997 






JP 


A2 


9215688 


19-08-1997 


WO Al 00030538 


02-06-2000 


EP 


Al 


1047338 


02-11-2000 



For more details about this annex see Official Journal of the European Patent Office, No. 12/82. 



16 



